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ABSTRACT 

The first stars ended the cosmic Dark Ages and created the first heavy elements necessary for the 
formation of planets and life. The properties of these stars remain uncertain, and it may be decades 
before individual Pop III stars are directly observed. Their masses, however, can be inferred from their 
supernova explosions, which may soon be found in both deep-field surveys by the James Webb Space 
Telescope (JWST) and in all-sky surveys by the Wide Field Infrared Survey Telescope (WFIRST). 
We have performed radiation hydrodynamical simulations of the near infrared signals of Pop III pair- 
instability supcrnovae in realistic circumstcllar environments with Lyman absorption by the neutral 
intcrgalactic medium. We find that JWST and WFIRST will detect these explosions out to z ~ 30 
and 20, respectively, unveiling the first generation of stars in the universe. 

Subject headings: early universe - galaxies: high-redshift - stars: early-type - supernovae: general - 
radiative transfer - hydrodynamics - shocks 



1. INTRODUCTION 

The cosmic Dark Ages ended with the formation of 
the first stars in 10 5 -10 6 M cosmological halos at 
z ~ 20-30. Besides ionizing the universe with in- 
tense UV flux and forming the first heavy elements, 
Pop III stars also determined the luminosities and spec- 
tra of primeval galaxies and may have been the origin 
of the supermassive black holes (SMBHs) at the cen- 
ters of most massive galaxies today. Unfortunately, be- 
cause they lie at the edge of the observable universe, 
individual Pop III stars will likely remain beyond the 
reach of direct observation for d ecades to com e de- 
spite their enormous luminosities (jSchaererl I2002D and 
new near infrared (N IR) observatories such as JWST 
(jGardner et all I2006Q and the Thirty Meter Telescope 
(TMT). Past numerical simulations have suggested that 
Pop III star s form in isolation in halos, with masses of 
30-300 Mrr. (iNakamura fc Umemuiall2001l iBromm et al.1 



[250l lAbel et all 120021: IP 'Shea fc Normanl I2007IL but 
mor e recent calculatio ns point to the possibility of bina- 
ries (|Turk et al.ll2009() . an d perhaps small swarms of less 
mass i ve 20-40 stars dStacv etall 120101: fClark et ak 



2011 



IGreif et al.ll20ll IHosokawa et al.ll201lHOreif et al 



20121 ). However, in spite of their increasing sophistication 



numerical simulations cannot yet constrain the propcr- 
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ties of the first stars, as they do not realistically bridge 
the gap between the formation and fragmentation of pro- 
tostellar disks and their destruction by nascent stars up 
to a Myr later. 

It may soon be possible to deduce the properties of 
Pop III stars from their supernova (SN) explosions. Pri- 
mordial stars from 10-40 M© are thought to die in 
core-collapse SNe while 140-260 M Q stars explode as 
pair-instability (PI) SNe, extremely energetic thermonu- 
clear explosions with up to 100 times the energy of 
Type la or II SNe rfHeger fc Wooslevl[200l IWhalen et al l 



2008bl Uoggerst et all 120101: iHeeer fc Wooslevl l201oT 
Joggerst fc Whalenll2011| ). Pop III stars above 100 M 



encounter the pair instability after central carbon burn- 
ing, when thermal energy creates electron-positron pairs 
rather than providi ng pressure support against collapse 
(|Barkat et al.lll967r i. Their cores subsequently contract, 
triggering explosive thermonuclear burning of G and Si. 
Above 140 M Q the energy released completely unbinds 
the star, with no compact remnant. At 260 M Q the 
core becomes so hot that alpha particles are photodisin- 
tcgrated into free nucleons, consuming as much energy 
per unit mass as was released by all preceding burn- 
ing, and the star collapses instead of exploding. PI SNe 
such as SN 2 007bi have now been discovered in the lo- 
cal universe (|Gal-Yam et al.l l2009f ). and the recent dis- 
covery that rotating Pop III stars can die as PI SNe at 
masses as low as 65 Mq may increase the number of 
such explosions at high redshif ts by a factor of a few 
(IChatzopoulos fc Wheeled [20l2h . 

Pop III PI SNe could be ideal probes of the primeval 
universe because they are 100,000 times brighter than 
their progenitors and the early galaxies in which they 
reside, and might be found by JWST or in all-sky sur- 
veys by WFIRST or the Wide-field Imaging Surveyor 
for High-Redshift (WISH). However, to determine if 
these events can be seen from Earth, fully radiation- 
hydrodynamical simulations of PI SNe must be per- 
formed in realistic circumstcllar envelopes, whose opac- 
ity and interaction with the shock crucially shape the 
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TABLE 1 

Pop III PI SN PROPERTIES (masses ARE IN Mq) 



M* 


M He 


R (10 14 cm) 


E (10 51 erg) 


M 6 6 Ni 


150 


72 


1.62 


9.0 


0.07 


175 


84.4 


1.74 


21.3 


0.70 


200 


96.7 


1.84 


33 


5.09 


225 


103.5 


3.33 


46.7 


16.5 


250 


124 


2.25 


69.2 


37.9 



temperatures and spectra in the source frame, and thus 
the flux that is eventually redshifted into the NIR. Ly- 
man absorption by neutral hydrogen prior to the era of 
rcionization, which can absorb or scatter light from these 
ancient supernovae out of our line of sight, must also be 
taken into account for explosions at the earliest epochs. 
Previo us simulations of PI SNe have not included these 
effects (jScannapieco et all 120051 iKasen eTaI]|20lll ), and 
estimates of the redshift at which they can be detected 
have either been confined to z < 15 (the era of primi- 
tive galaxies rather than the epoch of first light) or have 
been atte mpted at z ~ 30 with only very approx imate 
methods (pan et al.ll2012al lb1; iHummel et al.ll2012[ ). 

We have performed numerical simulations of the NIR 
light curves of Pop III PI SNe. In § 2 we describe our nu- 
merical method and in § 3 we present source frame light 
curves and spectra for 5 PI explosions. In § 4 NIR light 
curves for these events are shown for several redshifts, 
and we discuss their detection limits in redshift. 

2. NUMERICAL METHOD 

Our calculations are done in four stages. First, we 
model the evolution and initial explosion of 150, 175, 
200, 225, and 250 M R Pop III stars in th e Kepler code 
(| Weaver et all 1 19781 iWooslev et all 120021 ) . Second, the 
shock is propagated through the interior of the star, its 
surface, and then out into the surroundin g medium with 
the ra diation hydrodynamics code RAGE (iGittings et all 
120081 ). We model radiation flows with grey flux- limited 
diffusion and evolve matter and radiation temperatures 
separately because they are not always in equilibrium. 
This phy sics enables us to better c apture shock break- 
out (e.g. IMatzner fc McKeel fl999t iNakar fc SarH l27)ToT : 
IPiro et al boiOHKatz et al as discussed in greater 

detail in (jFrev et al.ll2012[ )7 Energy deposition in the gas 
due to the radioactive decay of 56 Ni, which powers the 
light curve at later times, is also included. 

We surround the star with a power-law wind profile 
with a total mass of 0.1 M , in keeping with the general 
consensus that massive zero-metallicity stars do not ex- 
hibit much mass loss due to their pristine atmospheres 
(jVink et all 1200 ll ). The wind profile is smoothly joined 
to the H II region density profile of a massive Pop III star 
calculated with the ZEUS -MP code (jWhalen et alJl2004t 
iWhalen fc Normanl 120061) . ZEUS-MP also captures the 
ionization state of the wind due to UV radiation from the 
star, and therefore its opacity to radiation from the SN. 
Because all five stars become red hypergiants and cease 
to emit ionizing radiation in the final few hundred kyr of 
their lives, the wind completely recombines by the time 
they explode while the low-density H II region beyond it 
remains ionized. 

Third, we construct spectra fro m the RAGE bla st pro- 
files with the SPECTRUM code (jFrev et al.ll2012l ) using 
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Fig. 1. — Source frame bolometric luminosities for the 5 PI SNe. 
The general trend of higher luminosity with progenitor mass is due 
to the larger explosion energies and mass of 5 °Ni that is synthe- 
sized. 

atomic opacities fr om the Los Alamos OPLIB database 
(jMagee et al.lll995D . SPECTRUM calculates continuum 
and line luminosities, absorption lines imprinted on the 
flux by the ejecta and envelope, and Doppler shifts and 
time dilation due to relativistic expansion. Finally, we 
convolv e our spectra with Lyman ab sorption by neutral 
clouds (|Madaulll99l ISu et all 120 111) , cosmological red- 
shifting, and filter response functions to calculate NIR 
light curves. We summarize the properties of our explo- 
sions in Table 1. 

3. SOURCE FRAME LIGHT CURVES / SPECTRA 

Bolometric source frame luminosities are shown for all 
five explosions in Figure [1] The shock is not visible to 
an external observer before breaking through the surface 
of the star because electrons in the intervening layers 
scatter the photons. When it reaches the surface, the 
shock abruptly accelerates in the steep density gradient, 
becoming even hotter and releasing a brilliant pulse of 
photons with peak luminosities above 10 47 erg s -1 , or 
2500 times the luminosity of the Milky Way galaxy. The 
transient is mostly X-rays and hard UV, and its duration 
is related to the radius of the star, since photons emit- 
ted simultaneously from its poles and equator reach an 
observer at times that are separated by the time it takes 
light to cross the star. The temperature of the shock at 
breakout is ~ 10 6 K. After the initial transient subsides, 
the luminosities gradually decline from 10 46 erg s _1 to 
10 42 erg s _1 over 3 yr as gamma rays from the decay of 
56 Ni formed during the explosion downscatter in energy 
and diffuse out through the ejecta. 

Pair-instability SN light curves are powered at early 
times by the conversion of kinetic energy into thermal en- 
ergy by the shock, and are thus far more luminous than 
Type la and II SNe because they have much higher ex- 
plosion energies. At later times they are primarily pow- 
ered by radioactive decay, and they remain much brighter 
than other SNe because they have synthesized more 56 Ni: 
up to 50 M compared to 0.4-0.8 M and < 0.3 M Q in 
Type la and II SNe, respectively. They are luminous for 
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Fig. 2. — Spectral evolution of the 225 Mq pair instability su- 
pernova. The fireball spectra are for 2.91 X 10 5 s (3.4 days, black), 
6.13 X 10 5 s (7 days, blue), 9.74 X 10 6 s (113 days, red) and 9.70 
X 10 7 s (1123 days, green). 

longer times (3 yr instead of 3-6 mo for Type la and II 
SNe) because radiation diffusion timcscales in their more 
massive ejecta are longer: 



(1) 



Here, k is the average opacity of the ejecta, M e j is the 
mass of the ejecta, and E is the explosion energy. All 
five light curves slightly rebrighten at 4 months when 
the photosphere (the surface from which photons stream 
freely into space) recedes into the hot 56 Ni layer deep in 
the ejecta. This emission peaks at ~ 3000 A, and its 
luminosity increases with the mass of 56 Ni, and hence 
the mass of the progenitor. Rcbrightcning occurs later 
in more massive stars because more time is required for 
the photosphere to sink into the 56 Ni in the frame of the 
shock. 

Figure [2] shows spectra for the 225 M Q explosion at 
four times ranging from shock breakout to 3 yr. As the 
fireball expands it cools, and its spectral peak advances 
to longer wavelengths. The breakout transient ionizes 
the envelope shrouding the star, as seen at 2.9 x 10 5 s 
(3.4 days) when the spectrum has far fewer absorption 
lines than at later times when the shock cools and can 
no longer fully ionize the wind. As the envelope recom- 
bines it imprints many more absorption lines on the spec- 
trum, especially at shorter wavelengths where line blan- 
keting effectively shears off the spectrum. Consequently, 
the common practice of fitting blackbodics to bolomctric 
light curves to approximate shock spectra can seriously 
overestimate the NIR flux in the observer frame, because 
such profiles have a significant amount of flux at short 
wavelengths that is actually removed by the envelope and 
ejecta. On the other hand, our bolometric luminosities 
are an order of magnitude higher at shock breakout than 
previous PI SN light curves because the shock crashes 
into the wind and heats it, and because radiation and 
matter are not forced to be in thermal equilibrium in 
our calculations. By 3 yr most of the photons in the 



24- 



26 



o 28 



E 30 




< 



32 



1148 A, no Ly at 
1432 A 
34 - 1148A,Lyabs 



36 . 



200 



400 600 
Time (days) 



800 



1000 



Fig. 3. — NIR light curves with and without Lyman absorption at 
z = 30. Blue: 3.56 fim (1148 A in the source frame), no absorption. 
Green: 3.56 /im with absorption. Red: 4.44 (im (1432 A in the 
source frame). 

spectrum are below the ionization limit of hydrogen and 
pass through the envelope, although resonant scattering 
continues to imprint absorption lines on flux below this 
limit. 

4. JWST AND WFIRST DETECTION LIMITS 

At high redshift, Lyman absorption will absorb most 
luminosity at wavelengths far blueward of 1216 A. How- 
ever, it could be possible that the fireball is brighter a 
little blueward of the Lyman limit rather than a little 
redward, even with absorption. To determine in which 
JWST NIRCam filter an explosion is brightest at a given 
redshift, we calculate its NIR signal at wavelengths above 
and below the Lyman limit in the observer frame. At z = 
30 we find that although all five SNe are brighter just be- 
low 1216 A in the source frame, IGM absorption reduces 
these luminosities to magnitudes below those just red- 
ward of the limit for which there is no absorption, as we 
show for the 225 M PI SN in Figure El The explosion 
is brighter at 1148 A in the source frame (3.56 fxm in 
the observer frame) at early times than at 1432 A (4.44 
/jm in the observer frame), but after IGM absorption is 
taken into account the event is brighter at 4.44 /zm in 
the observer frame. 

In the left panel of Figure 0] we show NIR light curves 
for all five Pop III SNe at 4.44 /an in JWST at z = 30, 
or 100 Myr after the big bang. The detection limit for 
this filter is AB magnitude 32, so all five explosions will 
be visible for over a thousand days. Their light curves do 
not cross above the detection threshold until 50 days (or 
1.5 days in the source frame at this redshift), so the most 
luminous component of the explosion, shock breakout, 
cannot be seen from Earth. The X-rays and extreme UV 
in the pulse are all redshifted into the UV and absorbed 
by neutral hydrogen. 

The NIR flux in the observer frame exhibits signifi- 
cantly more variability than the bolometric flux because 
of the expansion and cooling of the fireball, and this is 
essential to identifying these events as explosive tran- 
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Fig. 4. — NIR light curves for five Pop III PI SNe. Left: NIR magnitudes at 4.44 fim at z = 30, 100 Myr after the big bang. The horizontal 
dashed lines at mag 32 and 27 are JWST and WFIRST detection limits, respectively. The explosions remain above the NIRCam detection 
limit for over 1,000 days. Right: NIR light curves at 2.0 (im at z = 15 (solid) and z = 20 (dotted). 



sicnts. The flux rises much more quickly than it falls, 
so it is easiest to detect Pop III PI SNe in their ini- 
tial stages. Nonetheless, they exhibit sufficient variabil- 
ity over likely protogalactic survey times of 1-5 yr to 
be identified at later times as well. This variation over 
timescales of ^1000 days will also allow them to be dis- 
criminated from primeval galaxies, in spite of their prob- 
able overlap in color-color space. Although recent cal- 
culations show that roughly one Pop III PI SN will be 
present in any given J WST deep field dTrenti et alj|2009t 
iHummel etaLl 120121 : I Johnson etaD l2012t ). far greater 
numbers could be found in proposed all-sky NIR sur- 
veys with WFIRST. With a sensitivity of AB magni- 
tude 27 at 2.2 /im, WFIRST will see such explosions 
at z = 15-20, as we show in the right panel of Fig- 
ure |U Such redshifts may be optimal for detecting PI 
SNe because of the rise of UV background fields from 
the first generations of Pop III stars. This background 
is thought to delay primordial star formation in halos 
until they have grown to larger masses, forming more 
mass i ve stars at slightly lower redshifts (jWise k, Abell 
120071: lO'Shea fc Normanl I2008D . Star formation rates 
extracted from large-scale cosmological simulations of 
early structure forma tion that include UV backgrounds 
(| Johnson et all 120121 ) suggest that WFIRST will detect 
up to 10 3 PI SNe per year at 15 < z < 20. Detections of 
this volume will enable the first majo r surveys of the Pop 
III initial mass fun ction (IMF) (see lTanaka et al.l 120121 : 
iMoriva et al.l [201 2, for discussions on strategies for de- 
tecting superluminous Type Iln SNe in all-sky surveys). 
Our calculations definitively establish that JWST and 



WFIRST will be able to detect Pop III supernovae at 
the era of first light. In addition to unveiling the nature 
of primordial stars and constraining scenarios for early 
cosmological rcionization and chemical enrichment, de- 
tections of Pop III supernovae at slightly lower redshifts 
{z ~ 10-15) will probe the era of primitive galaxy for- 
mation, revealing their positions on the sky when they 
might otherwise have eluded detection by JWST. The 
discovery of the first cosmic explosions will be one of 
the most spectacular results in extragalactic astronomy 
in the coming decade, and open the first window on the 
end of the cosmic Dark Ages. 
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